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Abstract 

A linear relationship between the Hubble expansion parameter and the time derivative 
of the scalar field is assumed in order to derive exact analytic cosmological solutions to 
Einstein's gravity with two fluids: a barotropic perfect fluid of ordinary matter, together 
with a self-interacting scalar field fluid accounting for the dark energy in the universe. A 
priori assumptions about the functional form of the self-interaction potential or about the 
scale factor behavior are not neccessary. These are obtained as outputs of the assumed linear 
relationship between the Hubble expansion parameter and the time derivative of the scalar 
field. As a consequence only a class of exponential potentials and their combinations can 
be treated. The relevance of the solutions found for the description of the cosmic evolution 
are discussed in some detail. The possibility to have superaccelerated expansion within the 
context of normal quintessence models is also discussed. 



Dark energy or missing energy is one of the contemporary issues the physics community 
is more interested in due, mainly, to a relatively recent (revolutionary) discovery that our 
present universe is in a stage of accelerated expansion Jl|], that was preceded by an early 
period of decelerated expansion^. This missing component of the material content of the 
universe is the responsible for the current stage of accelerated expansion and accounts for 
2/3 of the total energy content of the universe, determining its destiny 0. This is a new 
form of energy with repulsive gravity and possible implications for quantum theory and 
supersymmetry breaking 0. 

A self-interacting, slowly varying scalar field, most often called quintessence, has been 
meant to account for the dark energy component. In any successful model of this class, the 
scalar field energy density should be subdominant at high redshift (in the past) and domi- 
nant at low redshift (at present and in the future) in order to agree with observations |fj. [3j. 
A variety of self-interaction potentials for the quintessence field has been studied. Among 
them, the simplest exponential potential (a single exponential) model is unacceptable be- 
cause it can not produce the transition from subdominant to dominant energy density 0. 
Other combinations of exponential potentials have been also studied||, 151 EL H, ]8[]. Com- 
binations of exponentials are interesting alternatives since these arise in more fundamental 
(particle) contexts: supergravity and superstring|J, where these types of potentials appear 
after dimensional reduction. 

In most cases the occurrence of a self-interaction potential for the scalar field makes diffi- 
cult to solve analytically the field equations, although some techniques for deriving solutions 
have been developed. In Ref.|T(|, for instance, the form of the scale factor is given a priori 



and, consequently, the self-interaction potential can be found that obeys the field equations. 
This method has been repeatedly used[ll|, [12|] as well as interesting applications of it [|13| . 



However there are cases when exact solutions can be found once the form of the potential 
is given ||], In other cases some suitable relationship between the self-interaction poten- 
tial and the scalar field kinetic energy is assumed ||14|| . The assumption that the scalar field 



energy density scales as an exact power of the scale factor has also been used for deriving 
solutions |15| . 



In this paper we explore a linear relationship between the Hubble expansion parameter 
and the time derivative of the scalar field to derive exact analytic cosmological solutions to 
Einstein's gravity with two fluids: a barotropic perfect fluid of ordinary matter, together 
with a self-interacting scalar field fluid accounting for the dark energy in the universe. A 



similar (but more general) relationship oc H has been used in Ref. Jl6j to find general 
solutions in N + 1-dimensional theories of gravity with a self-interacting scalar field and 
also to derive general solutions when we have a self-interacting scalar field with exponential 
potential plus a free scalar field.0 The assumed relationship between the Hubble parameter 
and the time derivative of the scalar field is suggested by an implicit symmetry of the field 
equations. We will be concerned with flat Friedmann- Robert son- Walker (FRW) cosmologies 
with the line element given by: 



*A method where a linear relationship between the field variables and/or their derivatives is assumed has 
been also used in p7[ to derive 4d Poincare invariant solutions in thick brane contexts. 
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ds 2 = -dt 2 + a(t) 2 S ik dx l dx k , (1) 

where the indexes i, k — 1, 2, 3 and a(t) is the scale factor. 

We point out that it is not neccessary to make any a priori assumptions about the 
functional form of the self-interaction potential or about the scale factor behavior. These 
are obtained as outputs of the assumed linear relationship between the Hubble expansion 
parameter and the time derivative of the scalar field, once one integrates the field equations 
explicitely. As a consequence only a class of exponential potentials and their combinations 
can be treated. However, this is not a serious drawback of the method since, as pointed out 
above, exponential potentials are of prime importance in dark energy contexts accounted for 
by a quintessence field. The relevance of the solutions found for the description of the cosmic 
evolution will be discussed in some detail. We use the system of units in which 8nG = c = 1. 

The field equations are: 

3H 2 = Pm + l -4> 2 + V, (2) 
2H + 3H 2 = (1 - 7 )p m - + V, (3) 

4> + 3H0 = -V, (4) 

where 7 is the barotropic index of the fluid of ordinary matter and H = ^ is the Hubble 
expansion parameter. The dot accounts for derivative in respect to the cosmic time t mean- 
while the comma denotes derivative in respect to the scalar field 0. The energy density of the 
ordinary matter (cold dark matter plus baryons and/or radiation) is related with the scale 
factor through p m = po/y^ -37 , where po,7 is an integration constant coming from integrating 
the conservation equation. Let us combine equations (^) and (|3j) to obtain 

H + 3H 2 = <1 ~^ L p m + 2V. (5) 

An implicit symmetry of the left hand side (LHS) of equations (|^) and (j^)is evident under 
the change H — > k(p. I. e., if one assumes a linear relationship between the Hubble expansion 
parameter and the time derivative of the scalar field; 

H = k(j), =^> a = e**, (6) 

where k is a constant parameter, the LHS of equations (f|) and (|]) coincide up to the factor 
k. We obtain a differential equation for determining the functional form of the potential V: 

V' + \V= ^Po, 7 e- 3 ^. (7) 

Explicit integration yields the following potential which is a combination of exponentials: 
V = £oe _I<?i + 6fc 2 2~l 2 Po>7 e ~ 3fc7 ^' where £0 is an integration constant. This potential can be 
given in terms of the scale factor if one considers equation (^): 
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V = toa'V + m^Po,^. (8) 

An interesting feature of this potential is that it depends on the type of ordinary fluid 
which fills the universe. Otherwise, it depends on the barotropic index 7 of the matter fluid. 
This fact implies some kind of interaction between the ordinary matter and the quintessence 



field much like the interacting quintessence studied in Ref . [£L8|] . This potential supports both 
an early inflation stage and a late time accelerated expansion through a stage of decelerated 
expansion. By substituting (§) back into Eq. (||), we can rewrite this last equation in the 
following form: (a/a) 2 = g^^ Po^a" 37 + Jp^-a -1 ^ 2 , or, if one introduces the following time 
variable dt = a l ^ 2k2 dr, 



where now the dot means derivative in respect to the new time variable r. This equation 

f a -1 da 

a 3 7 /2-l/2fc2_l da 



can be written in an integral form: / , a lrfa = r + To, or equivalently: 

V Aa 1 /fc 2 -37+_B 



y/A + Ba*y+V* 2 



r + r , (10) 



where A = gp^-^Avyj B = J^rj, and r is an integration constant. The integral in the LHS 
of Eq. ( |T(1| ) can be explicitely taken to yield the following expression for the scale factor as 
function of the time varaible r; 

a(r) = a {sinh[/i(r + ro )]} 2fe2 /(3^7-i) ; 



where a = [^^( 2 ^ I )] fe2/(3fc27 " 1) and fi = The Hubble expansion parameter 

can then be easily computed: 

Ik 2 

R2 ^ = ^37^ r )~ lA + *>)]■ (12) 

For the purpose of observational testing of the solutions it is useful to look for another 
magnitudes of astrophysical interest. Among them, the scalar field energy density: 



l -4> 2 + V = 3H 2 - p m , (13) 



the scalar field state parameter: 



\<p 2 + V 3k 2 (l-ll m ) 

where Q m = p m /3H 2 is the density parameter for ordinary matter that is related with the 
scalar field density paremeter through (just another writing for the field equation (Q)), 
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= 3^ = 1 - O m , (15) 
and last but not the less, the deceleration parameter q = —(1 + H/H 2 ): 

q = - 1 + W + Y Qm - (16) 

While deriving equations (|13|)-(|T6|) we have used the field equations (0), ([3]), (Q) and 
their combinations. Assuming the linear relationship (|j) between the Hubble parameter and 
the time derivative of the scalar field means that we are introducing a new free parameter, 
however, this parameter k can be assumed to be a known function of the other free parameters 
of the theory (r and £ ) or °f the barotropic index, etc. Another possibility, perhaps the 
most promising, is to choose a value for k according to the best fitting of the model to the 
observations. 

We now investigate the full space of solutions in k parameter without assuming apriori a 
definite form for it as function of the other parameters r and £o- 
A) k 2 < => k 2 = -k 2 > 

This case implies a negative kinetic energy for the scalar field since, due to @, 2 = 
-H 2 /k 2 . A kinetic energy term with the wrong sign has been advocated in Ref. [HJ to 
produce superaccelerated expansion. In this case the state parameter = — ^^Hyk 2 < — 1- 
I. e., in spite of contrary claims |20], |21] superaccelerated expansion and, correspondingly, 



superquintessence can be obtained within the frame of normal quintessence models like this. 
As it will be seen later on, from the observational point of view, for large k 2 (small negative 
kinetic energy contribution), superquintessence is undistinguishable from the case with usual 
accelerated expansion. We want to remark that, although the scalar field kinetic energy is 
negative, its energy density could be positive if 2V > H 2 /k 2 . Violation of the known 
energy conditions is obvious since = V$l P<t> < — 1 =^ P<t> + P<t> < US 



In terms of the positive parameter k 2 , the potential (H) takes the following form: 

V = £ a 1/P - -r^r^-^Po '■ (17) 
^ 2(3P7 + iy ' 7 V ; 

The scale factor behavior can now be rewritten in the following form: 

a(r) = a {sinh[/i(r + r )]} 2 ^ 3p ^, (18) 

2 7+l V 2£ 



where a = [^^(^ )] P/(3p7+1) and ji = J (3 g]+ 1 1)2 ^. For the Hubble expansion param- 



eter one has: H 2 {t) = 6 |2^ 1 ^ a(r) 1//fc2 coth 2 [/x(r + To)], meanwhile for the other magnitudes 
of observational interest: 

fi m = J 2 {cosh[^(r + r )]}- 2 . (19) 
For the state parameter (equation of state) and the deceleration parameter we have; 
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"* = - 1 - i»(i-a m y q = - {1 + W ) + T n "" (20) 

respectively. From this last pair of equations one sees that, for any finite k 2 > 0, the state 
parameter is more negative than —1, a phenomenon known as superquintessence (see 



Ref. [H], ^(| |2TJ] and references therein). Besides, an additional non vanishing negative 
contribution to the deceleration parameter is evident (compare the second equation in ( PU[ ) 
with (|16D), i. e., superacceleration is obtained. 



We shall notice that some undesirable properties of phantom energy |22|j are not present 
or are smoothed out in this case. For instance neither the scale factor of the universe nor 
the Hubble expansion parameter diverge in a finite (or infinite) amount of cosmic time 
(pointing into the future). Despite of the fact that, in general, the density of phantom 
energy p^ = 3k ^ffi 7 (sinh 2 \fl{r + r )] - ^)/(smh\fi(r + To )])6fe 2 7/(3fc 2 7+i) s i ow i y grows as 

the universe expands (for large k 2 » 1, p^ « £ ), we need an infinite amount of cosmic 
time to reach an infinite density state. 
B) k 2 > 

We should differentiate two ranges. For < k 2 < 1/37 the scale factor evolves periodi- 
cally: 

a(r) = a {sin[u(r + To)]} 2 ^ 3 ^" 1 ), (21) 

where a = [ ^-i) (^g^)] 2 ^ 3 ^' 1 ). However, in this case, for < 7 < 2, 3fc 2 7 - 1 < 
so one should have £0 < 0. This means that the self-interaction potential itself is negative, 
i. e., the scalar field energy density is negative. Besides q > always and this case is not of 
observational interest. 

For k 2 > I/37 the equations (|TTD, ( |T2"D are valid and the quantities of observational 
interest read; 



Mm = 3fc 3 p 7 1 { COsh [^( r + T o)]} 2 , 

^-- 1 + 3F(i 1 -n„,r « = -i + i/2* 2 + fa„. (22) 

Now we proceed to observationally test the solution found. We first compute the param- 
eters r , £o> an d k by the best fit to the observational evidence. These parameters can be 
fixed once one tries to fit the model with the experimental observations of SNla (we include 
the known bounds on f2o,m, go, etc.|2|, [|), including the recent observation of SN1997ff at 
redshift z = 1.7. We are left with a set in parameter space (t ,£o,&) of admissible values 
fulfilling with the observational bounds. From this set we choose a trio and then we compute 
the magnitudes of observational interest (equations (p!3|)-(^6|)). While doing the computa- 
tions we used the condition p ,m = n £,o- F° r illustration, in all figures in this paper, we 
chose the trio {n = 0.5, k = 10000, £ = 3) that agrees with the known bounds on fio,m, 
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flo,<f>, Qo, etc, and also we considered dust matter (7 = 1). For the case k 2 < 0{k 2 > 0) the 
present values of the state parameter and matter density parameter are uo^ = —1.00005 and 
Q,/, = 0.333341 respectively, meanwhile, for k 2 > 0, the present values of these parameters 
are cu = -0.99995 and fi = 0.333326. 

From fig.l, where the evolution of both Q m and & redshift z is shown, we see 
that the quintessence field 4> was subdominant in the past as required by nucleosynthesis 
constraints |24|, is comparable to the matter density parameter at present r = 0(z = 0) 
(according to present observations fio,m ~ 0.3 and Qq^ « 0.70) and will be dominant in 
the future. 

The dependence of the deceleration parameter on redshif is shown in fig. 2 for k 2 < 
0(k 2 > 0) and k 2 > 0(> I/37). The transition from decelerated expansion into accelerated 
expansion occurs approximately at z ~ 0.6 in both cases, in agreement with claims that the 
transition should be at z ~ 0.5 [[25|] and with recent observation of the SN1997ff at redshift 
z = 1.7, confirming a decelerated phase when the universe was a few seconds old. In fig. 3 
we show the distance modulus 5(z) as function of redshift both, computed according to the 
model (solid line for k 2 < and dashed line for k 2 > I/37) and the experimental curve 
(dots). Relative deviations are of about 0.7%. It is noticeable that both superquintessence 
and quintessence are undistinguisheable from the observational point of view, and for the 
values of the parameters considered here. 

The linear relationship explored in this paper allows, also, deriving solutions in Brans- 
Dicke and non-minimally coupled theories in general and this will be the subject of forth- 
coming papers. 

Summing up. We have derived exact analytic solutions to gravity theory minimally 
coupled to a self-interacting scalar field by assuming a linear relationship between the Hubble 
expansion parameter and the time derivative of the scalar field. This relationship is suggested 
by an implicit symmetry of the field equations. It induces a restriction upon the type of 
potentials one can deal with: a combination of exponentials that depends on the barotropic 
index of the matter fluid. We have derived exact analytic solutions and the entire range 
in k parameter has been explored. It has been shown that, contrary to common belief, for 
k 2 < superquintessence behavior is accomodated within this "normal" quintessence model. 
The solutions found have been observationally tested by using the distance modulus test 
(SNla observations where we included the known bounds on fi 0>m , <?o> e tc.|2], |3[]) and the 
SN1997ff observation at z = 1.7. The main features of present observational cosmology are 
satisfied within admissible accuracy levels for both k 2 < and k 2 > I/37, leaving a margin 
for superquintessence to effectively take place. 

We thank John D. Barrow and Luis P. Chimento for pointing out to us references Jl5| and 
T6| respectively. We acknowledge the MES of Cuba by financial support of this research. 
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Figure 1: The matter density parameter fl m (circled dotted line) and the scalar field density 
parameter (solid line) as functions of the redshift for 7 = 1, n — 0.5, k = 10000, and 
£0 = 3. We have used the condition p ,7 — n £,o- The parameters have been chosen such 
that, at present, ojq^ = —0.99995 and fi ,m = 0.333326. It is seen an early stage when 
the contribution from the quintessence field was subdominant. At present (z = 0) both 
contributions from dust and from the scalar field are of the same order (ilo,m = 1/3 while 
Qo,^ = 2/3). In the future the quintessence field will be dominant. 



Figure 2: Deceleration parameter q as function of the redshift z for 7 = 1 and values of 
the parameters n = 0.5, k 2 = 10000, and £0 — 3, that agree with the known bounds on 
f2 , m , Qo,<f>, Qo, etc. We have used the condition p ,7 — n £,o- The curves for k 2 < and 
k 2 > 0(> I/37) are shown toguether (we can not differentiate among them because they 
coincide). The transition from decelerated expansion into accelerated expansion occurs at 
z ~ 0.6 in agreement with known suggestions and with the observation of the SN1997ff at 
redshift z — 1.7 



Figure 3: Modulus distance vs redshift for 7 = 1 (n = 0.5, k 2 = 10000, £ = 3). The 
solid and the dashed lines represent the results of the theoretical model for k 2 < and 
k 2 > 0(> I/37) respectively. The dots account for the experimental data. A satisfactory 
agreement is achieved (the relative deviations are of approximately 0.7%). 
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